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Abstract 'Ghost Rods' are periodic structures in a two-dimensional flow that have 
an effect on material lines that is similar to real stirring rods. An example is a 
periodic island: material lines exterior to it must wrap around such an island, 
because determinism forbids them from crossing through it. Hence, islands act as 
topological obstacles to material lines, just like physical rods, and lower bounds on 
the rate of stretching of material lines can be deduced from the motion of islands 
and rods. Here, we show that unstable periodic orbits can also act as ghost rods, 
as long as material lines can 'fold' around the orbit, which requires the orbit to be 
parabolic. We investigate the factors that determine the effective size of ghost rods, 
that is, the magnitude of their impact on material lines. 
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1 Introduction 

Topological kinematics is the application of topology to chaotic advection in fluids. In 
two dimensions, braids are the natural mathematical construct to use for a topological 
analysis. Boylan d et al. used braids ve ry effectively to analyse the motion of stirring rods 
in vi scous flow ijBovland et al.l l2000(l and point vortices in ideal flow l|Bovland et al.l 
l2003(l . A braid is associated with the motion of the rods or vortices by plotting their 
trajectory in a space-time diagram: the resulting "spaghetti plot" is obviously a braid. 
Here, we shall not be too concerned with the precise mathematical properties of braids — 
the intuitive, capillary notion of what a braid resembles will suffice. 

Rods and points vortices share the common feature that they are topological obstacles 
to material lines in two dimensions. Of course, any fluid particle i s such an obstacl e, and 
recently one of us analysed braids formed by particle trajectories l)Thiffeaultl . l2005|) . The 
fact that particle orbits are topological o bstacles puts a lower bound on the topological en- 
tropy — the growth rate of material lines ijBovland et all I2OO0I: iNewhouse and Pignatarol 
Il99,'j^ . Imagine a material line that is initially linked around the topological obstacles 
under consideration (rods or fluid particles). Then as the position of these obstacles 
evolves the material line is dragged along, and as the particles braid around each other 
the material line must grow by at least a certain amount. The properties of the braid 
thus imply a lower bound on the growth rate of material lines in the fluid. 
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For time-periodic flows, the natural topological obstacles to study are fluid parcels 
associated with particular p eriodic orbits. Recently, we introduced to fluid mechanics 
the concept of "ghost rods" ijGouillart et al.li2005,) . We analysed the motion of material 
rods, stable periodic orbits (islands), and unstable periodic orbits from a topological 
perspective. We showed that periodic orbits associated with islands behave very similarly 
to material rods: they are large topological obstacles ploughing through the fluid, and 
material lines must get out of their way or else wrap around them. Periodic islands have 
the advantage of being easily identifiable visually, and can clearly be regarded as "rods." 
Figure ^ illustrates this: there is only one physical rod stirring the flow, but a ghostly 
second rod is clearly visible in the centre- left portion of the plot, aroun d which material 
lines are wrapped. Indeed, there is a regular island in that region l|Gouillart et al.[ 




Figure 1. A material line being stirred by a moving rod in a viscous fluid. The rod is the 
circle visible in the centre-right portion of the fluid, but observe that there is a rod-like 
structure in the centre-left portion. This is a periodic island that acts like a rod — a ghost 
rod. 




The foregoing descrip tion is topologic al in nat ure. The s ize of the rods is immaterial to 
the topological entropy (|Bovland et al.L 2000; F inn et al.L 1200,'^ . Of course, in practice 
their size matters a lot: if the rods are made smaller, so does the region of the fluid 
for which the topological entropy lower bound applies. In the limit of infinitely small 
rods, one might expect this region to shrink to zero. This is certainly true of physical 
rods: if they were made infinitely small, the fluid would not even notice their presence 
and nothing would happen, except in a vanishingly small region. There is currently no 
theory that gives the size of the affected region given the size of the rods and their path, 
but in practice it is observed (in viscous flows) that it is of the order of the size of the 
rods and the size of the region swept by their motion. 
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So much for physical rods. But what about ghost rods? As their name indicates, 
they have no material existence. However, since they behave much like physical rods, 
we may ask what is their effective size. That is, topologically a ghost rod is supposed 
to mimic a real physical rod, but how much of an impact does it effectively have on the 
surrounding fluid? For periodic islands, the answer is clear: the effective size of the ghost 
rod is the size of the island. Figure convincingly illustrates that, as far as material lines 
are concerned, there is a stirring rod of the size of the periodic island in the centre-left 
portion of the flow. 

For unstable periodic orbits, the answer is much less clear, since in principle ghost rods 
of this type have zero size. In this paper, we shall investigate the effective size of ghost 
rods associated with unstable periodic orbits. In fact, as we shall see, not all unstable 
periodic orbits can even be said to be ghost rods. Rather, only unstable periodic orbits 
of parabolic (as opposed to hyperbolic) type can hope to qualify as ghost rods. The local 
linear structure near an hyperbolic orbit prevents material lines from "wrapping" around 
the periodic point, so that it does not appear as a rod at all. For parabolic orbits of a 
certain type, the unstable manifold terminates at the periodic orbit, allowing material 
lines to wrap around the point without encountering the invariant manifold. Thus, the 
periodic orbit appears visually as a tiny rod, which is our criterion for considering periodic 
orbits to be ghost rods. 

2 Unstable Periodic Orbits 

In an incompressible flow, the linearised flow around an unstable periodic orbit can 
be one of two types. Figure |21 depicts the most common, called a hyperbolic orbit, 
or hyperbolic point if one is speaking of the Poincare section (stroboscopic map) of 
the time-periodic flow. There are two distinguished directions along which points 



respectively get away from or converge to the periodic orbit. These directions can be 
nonlinear ly extended into the unstable and stable manifolds of the periodic orbit. A 
sufficient condition for the periodic orbit to be hyperbolic is that its Floquet matrix have 
nondegenerate eigenvalues. (The Floquet matrix is obtained by linearising the system 
about the periodic orbit and integrating over a full period of the orbit, as we will do in 
SectionOlfor a specific flow.) Even though they are topological obstacles to material lines 
in the flow, such orbits can hardly be called ghost rods. This is because material lines 
must align with the unstable manifold of the periodic orbit, a phenomenon sometimes 




stable 



unstable 



Figure 2. The stable and unstable manifolds of a hyperbolic point. 
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referred to as asymptotic directionality l|Giona et all llQQ^lThiffeaulll 12004) . A material 
line cannot possibly fold around such a periodic orbit, since the unstable manifold goes 
straight through the orbit and appears linear in its neighbourhood. Hence, the periodic 
orbit does not "look" like a tiny little rod to the naked eye: it looks like any other 
point on the material line, and only a detailed knowledge of the velocity field allows its 
detection, usually by numerical means. We conclude that hyperbolic periodic orbits do 
not form "proper" ghost rods, since they cannot be detected visually. 

That leaves the second type of unstable periodic orbit: parabolic orbits. In that 
case the Floquet matrix has degenerate eigenvalue that must both be equal to unity, by 
incompressibility of the fluid. For parabolic points, we cannot deduce the behaviour of 
points near the periodic orbit by examining only the linearised system — nonlinear terms 
must be considered. As we shall see in the following section, a particular type of nonlinear 
structure gives rise to parabolic points that exhibit the appropriate behaviour for a ghost 
rod. 



3 Case Study: The Sine Flow 

We shall now illustrate the type of unstable periodic orbit that gives rise to ghost rod s 
by examining a specific system, namely the Zeldovich sine flow IjPierrehumber 
This is a nice system to work with because its Poincare map can be obtained analytically, 
and for special parameter values we can also determine some periodic orbits analytically. 
These orbits were exploited bv iFinn et alJ l|2n0,'i|) to show the presence of chaos. The 
sine flow is given by the velocity field 

I ,N J (0, sin27ra;), tit < t < (n + i)r; 

|(sin27ry, 0), (n + i)T < t < (n + l)r, 

where n is an integer. The equation x — u(x^ t) can then be integrated over one period 
to give the sine map 

Xn+i = x„ + i r sin27r2/„+i ; 

Vn+l = Vn + I T sin 2nXn , 

with Xn '■— x{nT). As an example, we will take t — 1, because then we can determine 
many periodic orbits analytically. For instance, there is a period-4 orbit starting at x^ = 
(1/12, 1/2), with iterates 

(1/12, 1/2) (7/12, 3/4) ^ (7/12, 1/2) ^ (1/12, 1/4) ^ (1/12, 1/2). (3.3) 



The initial location of this orbit is inside the small square in Fig. 3(a) which also shows 



a material line advected for a few periods of the sine flow. Figure 3(b) is a blow-up of 



the material line near this periodic orbit. Notice how the material line is sharply folded 



around the periodic orbit. In fact. Fig. 3(a) contains several such sharp folds. They 



are quite g eneric in chaotic flows , and are associa ted with regions of anomalously low 
stretching l)Liu and Muzzid. Il996t iThiffeaultl l2004|l . 
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(a) 



(b) 



Figure 3. (a) A material line stretched and folded by the sine flow. The parabolic 
periodic point at {x,y) = (1/12, 1/2) is shown boxed, (b) Blow-up of the periodic point. 
Note that the material line is folded near, but not quite tightly around, the parabolic 
fixed point. 

We are interested in the behaviour of the map (|3.2|l in the neighbourhood of this 
periodic orbit, so we define a variable X := x — Xq. Then, we can expand the map to 
second order in X , 



where the primes denote four iterations of the sine map, so that the periodic orbit has 
become a fixed point of the map (|3.4() at X = (0, 0). The periodic orbit is parabolic, as 
can easily be seen from the fact that the linear part of (|3.4(l (the Floquet matrix) has 
matrix representation 



which implies unit eigenvalues for the map. However, this matrix is not diagonalisable: it 
only has one eigenvector, (0 1)"^ (this can only occur for a parabolic point). We will see 
that it is this nondiagonalisable nature that allows the "folding" of material lines around 
the periodic point. In general, a matrix 3 has this property if {3 — I)'^ = 0, for 3 7^ I, 
which given that det3 = 1 is equivalent to tr3 = 2, with 3 7^ I. 



X' = X- 2nY + a^X^ + (3^Y^ XY, 

Y' ^Y +ayX^ + (3yY^ + Jy XY, 



(3.4) 




(3.5) 



After a linear transformation and a near-identity area-preserving quadratic transfor- 
mation, Eg. 13.41 can be brought into the form 

X' = X + Y + aXY , 

Y' = Y + a{\X^ + XY) , ^^'^^ 

where the coefficients are such that the map is area-preserving to linear order. (The 
transformation used to get to is not generally orientation-preserving.) As long as 
the linear part of the system is a Jordan block of the form (|3.5|) , we can transform the 
system to Eq. (I3.6|l . We have thus reduced the dynamics near the parabolic point to a 
one-parameter map (basically a Henon map), which we proceed to analyse. 



3.1 Invariant Manifolds and Dynamics Near the Origin 

We now want to find the shape of the unstable and stable manifolds of the fixed point 
of lESl) at the origin. Unlike hyperbolic fixed points, for a parabolic point the invariant 
manifolds cannot be determined solely from the linear part. Rather, we use the invariance 
property of the manifold: we parametrise the invariant manifold by {X,Yi^y{X)) and 
iterate the map, 

X' = X + Yin^{X) + aXY;n^{X), (3.7a) 
51„,(X') = Y,^,{X) + a {\X^ + XY,^,{X)) , (3.7b) 

where we wrote Y' — Y{t^^{X') since, by the invariance property, the iterated point must 
still belong to the invariant manifold. We can then substitute H3.7a|l into 13.7b(l . 

11„,(X + yinv(X) +aXY;^AX)) = >lnv(^) + a + XY,,„{X)) , (3.8) 

which is an equation that must be solved for Yi^„{X). We are interested in the small X 
form of the manifold, so we assume Y\^„{X) = aX^ and try to balance the leading order 
terms; 

a{X + aX^ + a aX^+^)^ = crX^ + a {^X^ + aX^+^) . (3.9) 

Where we go next depends on the magnitude of S. US — 1, we get the equation cr{l+(T) = 
a for the coefficients of the linear terms, which implies u = 0, an unacceptable state of 
affairs since then the quadratic term is unbalanced. If (5 < 1, we get the leading balance 

a{<7X^)^=aX^. (3.10) 

This can only be satisfied for 6 — 1, a contradiction, or 6 = 0, which again leads to an 
unbalanced quadratic term. Hence, our only choice is to take S > 1, which gives the 
leading-order balance 

(j{X + aX^)^ ^crX^ + ^aX^, (3.11) 

where we have kept an extra order, because the leading terms cancel after expanding the 
exponent, 

aX^ (l + aSX^-^) = aX^ + ^aX^, (3.12) 



6 



and we get finally 



25-1 



(3.13) 



yielding (5 = 3/2, ct = ±\Ja/i. Wc can thus write the shape of the invariant manifolds 
as 



to leading order, where we have written aX under the square root to show that X and a 
must have the same sign. Note an important fact: the manifolds exist only on one side of 
the X axis, in contrast to Fig. El where the manifolds must radiate from the fixed point 
in four directions. 

The two solutions for a correspond to the stable and unstable manifolds. We can 
determine which sign goes with which manifold by looking at an iterate of X in Eq. H3.7a|l . 



where we neglected higher-order terms (X^/^). Thus, for a > 0, which implies X > 0, 
the '+' solution takes X farther from the origin (unstable manifold), whilst the '— ' 
solutions takes the point closer to the origin (stable manifold). The situation is reversed 
for a < 0. 

Figure 01 shows a few iterates of horizontal lines under the action of the map (|3.6|l . 
The linear part of the map acts as a "shear flow" that sweeps the line around the origin, 
but the nonlinear terms prevent the line from crossing the unstable manifold. The net 
result is a line folded around the unstable manifold. This is why it is appropriate to 
refer to (|3.6|l as the 'folding normal form': inside every sharp fold of the flow lurks such 
a map.^ Since nearby material lines align with the unstable manifold of the periodic 
orbit, the folding is made possible by the one-sidedness of the unstable manifold: unlike 
hyperbohc points (Fig. EJ, the manifold does not traverse the parabolic periodic point, 
but instead terminates there. This allows material lines to wrap around the periodic 
point without encountering the invariant unstable manifold, which cannot be crossed. 

3.2 Curvature 

As time progresses the folds in the material lines in Fig.^jcome closer and closer to the 
periodic orbit. There seems to be no limit to how close they can come, which is consistent 
with the ghost rod having zero effective size. The best we can do is to characterise the 
effective size of the ghost rod by how quickly the curvature of the folds evolve. We shall 
now examine how the curvature of a material line evolves near the parabolic orbit. 

Consider a material line going through the origin of (ISU, as depicted in Fig.EI The 
tangent map of H3.6|l at the origin tells us how the tangent to the curve evolves. 




(3.14) 




(3.15) 




(3.16) 



^At the workshop, Stefano Cerbelli and Massimiliano Giona pointed out that their research 
also seems to support this. 
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Figure 4. Iteration of a few material lines by the map (|3.6() for a = 1. The lines fold 
around the unstable manifold (dashed curve). The X axis has been rescaled by 10~^, 
the Y axis by 10"^. 



where {5X SYY^ is the tangent. The second variation of H3.6() is 

(5^X'\ _ (\ 1\ (5'^X\ ( 5X5Y 
\5^Y')^\Q \)\5-^Y)^'^\(5Xf + 5X5Y 

For the case shown in Fig. 0] the initial tangent is parallel to (1 0)"^, which is an 
eigenvector of the matrix in H3.16|l : the tangent doesn't change. Given that 5^X and 5^Y 
are zero initially (the line is straight), we can solve H3.17|l for 5^Y , 

5'^Y = na{5Xf , (3.18) 



wher e n is the number of iterations 

K 



(3.17) 



The curvature of the line is given by ijLiu and Muzziol 
_ SXS^Y - SYS'^X 

~ Wsxf ■ 



8 



Now given the solution H3.18|l and the fact that 5Y = for all time, the curvature evolves 
as 

K. = na, (3.20) 

so the curvature of the material line grows linearly with time. This is verified by a 
calculation with the sine flow, for the periodic orbit H3.3|l : Figure O shows that the 
curvature of a material line anchored at the periodic orbit does indeed grow linearly 




5 10 15 



periods 

Figure 5. The evolution of the curvature of a material line passing through the periodic 
point shown in Fig.O with orbit given by Eq. (|3.3() . The curvature at the point increases 
linearly with the number of periods, showing that the line is getting more tightly folded 
around the parabolic point. 

with the number of iterations, and the slope is in perfect agreement with the results 
from the folding normal form. 

Note that this linear evolution of the curvature is not an artefact of our choice of 
orientation of the material line. If we choose the line to be orthogonal to the one in 
Fig. ^ we find that the tangent evolves as {5X SYY" = {n ^)^5Yq, which means that 
the tangent aligns with the direction of the unstable manifold for large n. The curvature 
will then grow asymptotically at the rate given by H3.2U|I . 

It is clear from Fig.^that the point of highest curvature is not at origin. Nevertheless, 
the increase in curvature is linear in the neighbourhood of the periodic orbit, and all 
material lines near the orbit eventually wrap around it, so the orbit can still be said to 
be acting as a ghost rod. 

We conclude that a measures the "size" of the rod: a higher a means that material 
lines converge towards the periodic orbit more rapidly, so that the ghost rod has a smaller 
apparent impact on the flow. Visually, a could be estimated by looking at the rate at 
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which material hnes "bunch up" near a periodic orbit, as in Fig. 0] but in practice this 
is quite difhcult. 

4 Discussion 

Of course, this paper is just a sketch of a theory for the size of ghost rods: a comprehensive 
theory remains to be developed. Rather, we tried to give an indication of the factors 
that influence a ghost rod's apparent size. 

The motivation behind this study, and the ghost rod framework in general, is to try 
to determine some stirring properties of a flow from visual cues. It is obvious that we can 
determine the size of physical rods by just looking at them. The effective size of ghost 
rods that are elliptic islands can also be determined visually, as is apparent in Fig. ^ 
Ghost rods associated with parabolic points are harder to identify: they typically occur 
inside sharp folds in material lines, as in Fig. Even if they are identified, measuring 
their effective impact on the flow is far from trivial: one can attempt to measure the 
evolution of curvature near the point, in the same spirit as in Section 13.21 or see how 
rapidly material lines bunch-up near the periodic orbit, effectively a measure of the 
coefficient a. 

The sharp folds observed in material lines are the spots where the stretc hing is weak- 
est, because there is us ually a competition between stretching and curvature l|Liu and Muzzi( 
ll99fil:lThiffeimiill2n04l) . Hence, folds are associated with inhomogeneities in the stretching 
field, and thus typically decrease the efficiency of stirring since uniformity is desirable. 
Knowing how fast the curvature grows (as measured by a) gives a hint of how much 
inhomogeneity a fold introduces, and thus quantifies its impact on the quality of mixing. 

The parabolic points may be the "relevant" ghost rods, i.e. the ones o n which one can 
const ruct a braid that captures exactly the topological entropy of the flow l|Goiiillart et al.L 
We have no proof of this assertion yet; however, since the folds determine the skele- 
ton around which a material line will wrap, these points certainly play a distinguished 
role. 
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